Evolution of the proton sd states in neutron-rich Ca isotopes 
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We analyze the evolution with increasing isospin asymmetry of the proton single-particle states 
2sl/2 and ld3/2 in Ca isotopes, using non-relativistic and relativistic mean field approaches. Both 
models give similar trends and it is shown that this evolution is sensitive to the neutron shell 
structure, the two states becoming more or less close depending on the neutron orbitals which are 
filled. In the regions where the states get closer some parametrizations predict an inversion between 
them. This inversion occurs near 48 Ca as well as very far from stability where the two states 
systematically cross each other if the drip line predicted in the model is located far enough. We 
study in detail the modification of the two single-particle energies by using the equivalent potential 
in the Schroedinger-like Skyrme-Hartree-Fock equations. The role played by central, kinetic and 
spin-orbit contributions is discussed. We finally show that the effect of a tensor component in the 
effective interaction considerably favors the inversion of the two proton states in 48 Ca. 
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I. INTRODUCTION 

Novel properties and new scenarios are expected for 
nuclei situated far from stability. The new generation 
of radioactive beam facilities will allow to answer many 
open questions about the pecularities of these unstable 
systems. One of the major issues in the physics of exotic 
nuclei is the study of shell structure and magicity evolu- 
tion when approaching the drip lines [l|, 0] . From a the- 
oretical point of view, two aspects have been underlined 
as mainly responsible for the evolution of single-particle 
energies far from stability, the one-body spin-orbit poten- 
tial which is strongly modified when the surface becomes 
more diffuse [l| and the tensor force between neutrons 
and protons in valence sub-shells Q . 

Recently, the TV = 28 shell closure has been experimen- 
tally analyzed in the 46 Ar(d,p) 47 Ar transfer reaction [J. 
A strong reduction of the neutron p spin-orbit splitting 
has been observed in 47 Ar with respect to the isotone 
49 Ca. Since p states are mainly localized in the interior 
of the nucleus, this strong reduction cannot be justified 
by the presence of a diffuse surface which would affect 
only high-Z states mainly concentrated at the surface. A 
theoretical analysis based on the relativistic mean field 
(RMF) approach has been proposed by Todd-Rutcl et 
al. @. It predicts a strong reduction of the spin-orbit 
splitting for neutron 2p-states in 46 Ar as compared to 
48 Ca. At Z = 20, the state 2sl/2 is usually located be- 
tween ld5/2 and ld3/2. In the RMF calculations of Ref. 
0, however, 2sl/2 is predicted less bound than ld3/2 in 
both 46 Ar and 48 Ca (2sl/2 - ld3/2 inversion). In this 
scenario, 2sl/2 is empty in 46 Ar and occupied in 48 Ca: 
thus, the proton density profile in 46 Ar presents a strong 
depletion in the interior of the nucleus. This reduction of 
the charge density in the center would be responsible for 
the modification of the spin-orbit in the nuclear interior 
and, hence, for the reduction of the neutron 2p-splitting. 



This problem of 2sl/2 — ld3/2 inversion of the proton 
states has been already analyzed by Campi and Sprung 
@ within the Hartree-Fock (HF) + BCS model with an 
interaction derived from a G-matrix 0. 36 Ar was found 
as a candidate for this inversion. Skyrme forces do not 
predict any inversion in this nucleus. It is thus worth- 
while to revisit the problem for other nuclei in this region 
of the nuclear chart in the framework of the Skyrmc-HF 
model. 

In this work, we analyze the evolution of the s — d pro- 
ton single-particle states in Ca isotopes and the possible 
2s 1/2 — ld3/2 inversions. We also present some com- 
parisons with the corresponding results obtained within 
RMF. We neglect pairing in our treatment since Ca iso- 
topes are proton closed-shell nuclei. We have checked 
that, within RMF the inclusion of neutron pairing does 
not modify in a significant way the evolution of the pro- 
ton states we are interested in. The only important effect 
due to pairing is the shift of the drip line towards heavier 
isotopes (for example, the drip line is shifted from 60 Ca 
to 76 Ca with the parametrization NL3 [8j). However, this 
aspect is not relevant in the present analysis which is not 
intended to make any prediction on the drip line position. 
We choose the Ca isotopes since experimental signals for 
the inversion phenomena have been found at least in one 
of these isotopes, 48 Ca: the ground state of 47 K (one 
proton less than 48 Ca) is 1/2+ with a large spectroscopic 
factor [§] and the single-particle spectrum of 48 Ca has 
been measured, the proton state ld3/2 being more bound 
than 2sl/2 by about 300 keV In our analysis, we 

explore all the contributions, kinetic, central, spin-orbit 
and tensor, which can modify the single-particle energies 
with increasing A and we show that not only the spin- 
orbit and tensor terms are determinant. The role of the 
central mean field term is in particular discussed. Within 
the models which predict the crossing between the two 
states, we show that this inversion occurs near 48 Ca as 
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well as in very neutron-rich nuclei close to the drip line. 

The article is organized as follows. In Sec. II we study 
the evolution with increasing A of the difference Ae be- 
tween the 2sl/2 and ldi/2 energies obtained within non- 
relativistic and relativistic approaches. In Sec. Ill wc 
concentrate on the non-relativistic case and perform a de- 
tailed analysis of the results. The different contributions 
to Ae arc isolated by analyzing them with the equivalent 
potential in the Schroedinger-likc HF equations. In Sec. 
IV the effect of the tensor force is estimated in the frame- 
work of the SLy5-HF [ll[ model. Finally, conclusions are 
drawn in Sec. V. 



II. EVOLUTION OF 251/2 AND 1D3/2 PROTON 
STATES WITHIN SKYRME-HF AND RMF 

We first perform a preliminary study with HF calcu- 
lations of 48 Ca using different Skyrme interactions. We 
then choose three representative forces: SkI5 [l2j which 
predicts a 2sl/2 — le?3/2 inversion with an energy differ- 
ence Ae of ~ 800 keV, SGII [13J which also reproduces the 
inversion (Ae ~ 200 keV) and SLy4 [ll[ for which there 
is no inversion. With the three selected parametrizations 
we have systematically analyzed the Ca isotopes from 
40 Ca up to the HF two-neutron drip line. We recall that, 
in the three considered Skyrme parametrizations there is 
no explicit tensor force. 

We show in Fig. I the difference Ae between the en- 
ergies of the proton states 2sl/2 and ld3/2 for the three 
Skyrme forces. The inversion takes place where Ae is 
positive. Corrections to the individual energies due to 
the coupling of single particle motion with collective vi- 
brations, which arc neglected in our treatment, should 
be expected (see, for instance, Ref. [HI)- However, by 
considering the energy difference instead of the individ- 
ual single-particle energies the effects of these corrections 
should be reduced, the coupling to vibrations having the 
tendency of shifting upwards the energies of both occu- 
pied states. 

We observe that the SLy4-HF calculations never pre- 
dict the 2sl/2-ld3/2 inversion. On the other hand, both 
SkI5-HF and SGII-HF predict this inversion around 48 Ca 
as well as for more neutron-rich isotopes starting from 
58 Ca up to the drip line. The HF two-neutron drip line 
is located at 82 Ca, 78 Ca and 60 Ca with SkI5, SGII and 
SLy4, respectively. The two experimental points for 40 Ca 
and 48 Ca are also included in the figure. The three sets of 
results globally present the same behavior. Indeed, in all 
three cases the quantity Ae starts from a negative value 
and increases from A=40 to A=48. This generates a 
2sl/2-ld3/2 inversion with SkI5 (in 46 Ca, 48 Ca and 50 Ca) 
and with SGII (in 48 Ca). Going from 48 Ca up to 52 Ca 
the states cross again with SkI5 and SGII, while the dis- 
tance between them increases with SLy4. Beyond 52 Ca 
Ae increases again with the three parametrizations. This 
generates another inversion with SkI5 and SGII starting 
from 58 Ca. We notice also that the nuclei for which Ae 
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FIG. 1: Difference between the energies of the 2sl/2 and 
ld3/2 proton states calculated with the Skyrme interactions 
SkI5, SGII and SLy4 for Ca isotopes. 



presents maxima or minima are the same for the three 
Skyrme forces. 

A natural question to ask is whether the above general 
trends are specific of the Skyrme-HF approach, ft is well 
known that the RMF approach gives a spin-orbit poten- 
tial whose (N-Z) dependence is somewhat different from 
that of Skyrme-HF models^]. We have performed RMF 
calculations with different parametrizations for the same 
set of Ca isotopes using the parametrizations DDME1 
[111, NL3 [I] and NLB2 [T3| - The latter one is chosen 
as an example of RMF model which does not lead to a 
2sl/2-ld3/2 level inversion. The calculated values of Ae 
are shown in Fig. 2 up to 60 Ca which is the two- neutron 
drip line isotope predicted by DDME1 and NL3. Glob- 
ally, we observe for Ae the same trend as that obtained 
within the non relativistic HF, with maxima and minima 
corresponding to the same nuclei, 48 Ca and 52 Ca. Since 
comparable trends are obtained in both non-relativistic 
and relativistic approaches we conclude that the calcu- 
lated evolution of 2sl/2 and ld3/2 states is a generic 
behavior. We can thus explore more in detail the results 
by considering only the non-relativistic case. 



III. ANALYSIS OF THE CONTRIBUTIONS TO 

Ae 



We now concentrate on the maxima and minima of 
Ae. They correspond to nuclei with neutron closed shells 
or sub-shells: the maximum at 48 Ca corresponds to the 
closure of the neutron 1/7/2 orbital whereas the mini- 
mum at 52 Ca corresponds to the filling of the neutron 
2p3/2 state. In the non relativistic Skyrme-HF model 
the radial HF equations can be expressed in terms of an 



3 




7 

6 
5 
4 

> 3 
u 

S 2 
EP 1 

o 

w 
-1 

-2 
-3 

-4 



. 1 1 1 ' i 1 

^ — \ 
- w" \ 


i i i i 
__ — — — ' " 


l 1 l 

__ — ■ _ 

- 




• • Kinetic 






■ ■ Central 






* ■ ■ ■ * Spin-orbit 













40 44 



48 



52 56 
A 



60 64 



72 



FIG. 2: Difference between the energies of the 2sl/2 and 
ld3/2 proton states in Ca isotopes calculated in RMF with 
the parametrizations DDME1, NL3 and NLB2. 



FIG. 3: Kinetic, central and spin-orbit contributions of Eq. 
© for SkI5 in 40 Ca, 48 Ca, 52 Ca and 70 Ca. 



energy-dependent equivalent potential V l J q : 



2m 
where 



dr 2 



ijj(r) 



V c 



1(1 + 1) 



ijj(r) 



7«(r,e)^(r) 



+ ■ 



-U l s i(r)- 



m*(r) 



(1) 
e, (2) 



with U%(r) = U so (r) x \j(j + 1) -1(1 + 1) - 3/4]. U so (r) 
is the spin-orbit HF potential and V^ ntr ' is 
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where Uo(r) is the central HF potentials and m*(r) is the 
effective mass [ll[ • For protons Uq includes the Coulomb 
potential. Up to a normalization factor the HF radial 
wave function <f> of energy e is related to the solution ijj 
of Eq.((l|) by the relation ?/> = (m* /m) 1 / 2 ^. From Eqs. 
CO)-© we can write Ae as 
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where T is the kinetic contribution. The 3 terms of the 
r.h.s. of Eq. Q - kinetic, central and spin-orbit - are 
plotted in Fig. 3 for the force SkI5 and the nuclei 40 Ca, 



48 Ca, 52 Ca and 70 Ca. Similar results are obtained with 
SLy4 and SGII. We mention that the mean value of the 
effective mass in the denominators of Eq. ([4]) has very 
little A-dependence from 40 Ca to 70 Ca. From Fig. 3, one 
notices that the spin-orbit and kinetic terms present a 
regular behavior as a function of A. Both of them are 
weakened with increasing isospin asymmetry favouring 
the inversion in very neutron-rich isotopes. The spin- 
orbit term is weakened because the neutron surface be- 
comes more diffuse with increasing A. In general, the 
kinetic energy of an orbital depends on the mean dis- 
tance between its single-particle energy and the bottom 
of the potential in the region where the wave function 
is localized. For the two states 2sl/2 and ld3/2 we can 
look at the difference eij — V^'(ro), where ro is the root 
mean square (r.m.s.) radius of the corresponding wave 
function. The evolution of V^ q with increasing A is gov- 
erned by two effects: i) the lowering of the proton po- 
tential due to the symmetry term; ii) the formation of a 
neutron skin which modifies the proton distribution by 
pulling it towards larger radii. The intensity of these two 
effects depends on the quantum numbers of the neutron 
orbitals which are filled and of the proton wave function 
under study. As an illustration, we consider 52 Ca and 
70 Ca. The r.m.s. radii ro and the values ey — VM(ro) are 
shown in Table I for the 2sl/2 and ld3/2 states. From 
52 Ca to 70 Ca the difference ey — V^Uro) is reduced more 
for ld3/2 (4.9%) than for 2sl/2 (0.4%). This analysis is 
confirmed by the evolution of the two r.m.s. radii. It is 
evident that, going from 52 Ca to 70 Ca, the ld3/2 wave 
function is more affected than 2sl/2 by the enlarging of 
the potential due to the formation of a thick neutron skin. 
This explains the increase of the kinetic contribution to 
Ae with the neutron number. We consider now the cen- 
tral term of Eq. Q which is responsible for the maxima 
and minima of Ae, and concentrate on the maximum at 
48 Ca. We mention that the major role played by the 
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A 


ro (fm) 
2sl/2 


ro (fm) 
l<23/2 


e-V%(ro) (MeV) 
2sl/2 


e-V%(r ) (MeV) 
ld3/2 


52 


3.71 


3.70 


23.11 


19.73 


70 


3.81 


3.89 


23.01 


18.76 



TABLE I: The values of ro and e — V^(ro), for the states 
2sl/2 and ld3/2 in 52 Ca and 70 Ca. The interaction is SkI5. 
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FIG. 4: Vo, Vi and V 2 calculated with SkI5 in 40 Ca, 48 Ca and 
52 Ca. 
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central term in modifying the single particle energies has 
been also underlined by Gaudefroy et al. [J]. We intro- 
duce the quantities Vb, V\ and V2 which correspond to the 
contributions of the 3 terms of Eq. (|3|). They are plotted 
in Fig. 4 for 40 Ca, 48 Ca and 52 Ca. It turns out that the 
term mainly affected by the neutron shell structure is Vo 
which contains the Hartrce-Fock potential. 

We can separate the energy contributions of the 
N=Z=20 core from those of the excess neutrons. For 
instance, the total nucleon density p is a sum of p cor e 
and Pexcess, and similarly for the other types of densi- 
ties. Then, for any HF quantity the core contribution is 
obtained by replacing in its expression the total densities 
by core densities whereas the neutron excess contribution 
corresponds to the rest. We show in Fig. 5 the core and 
neutron excess contributions to Vb- It is clear that the 
change of slope at 48 Ca is mainly due to the neutron ex- 
cess contribution. We have further verified that the term 
mainly responsible is the £0 term of the Skyrme force. 
The density-dependent term (t3 term) is also sensitive to 
the neutron shell structure but with an opposite behav- 
ior reducing the effect due to the to term. Hence, the 
main parameters which influence the behavior of Ae are 
Xq and to as well as £3, £3 and a. We have checked that 
the role played by the other terms of the HF potential is 
negligible. 

To complete our analysis we consider separately the 
two single-particle energies 2sl/2 and ld3/2. We have 



> 



b 3^ 

O 

W 2- 



( 3 



1 1 1 

m——— ~~~~~~ 


1 1 1 1 1 ' 


1 

- 


- 












- 


* ■ ■ ■ * n excess contrib. 


- 




■ — ■ core contrib. 






...*.. 




J, 








40 



42 



44 



46 
A 



48 



50 



52 



54 



FIG. 5: Core and neutron excess contributions to Vo for SkI5 
in 40 Ca, 48 Ca and 52 Ca. 



verified that the maximum of Ae at Ca is mostly due to 
the energy of the 1g?3/2 state which decreases less rapidly 
from 48 Ca to 52 Ca than from 40 Ca to 48 Ca. This behav- 
ior is explained in Fig. 6. In the top panels the neutron 
excess contribution Phf to the HF potential due to the 
to and £3 terms is plotted for 44 Ca, 48 Ca and 52 Ca (see 
e -g- [fl| for the expressions in terms of the Skyrme force 
parameters). In the middle panels we show the squares 
of the ld3/2 and 2sl/2 radial wave functions multiplied 
by r 2 . In the bottom panels the products Phf\<^"t 2 
are shown. Since Phf is negative (the to contribution 
is negative and the largest in absolute value) the proton 
single-particle energies are lowered with increasing N-Z, 
as expected. We observe that, in 44 Ca and 48 Ca the po- 
tential related to the neutron excess (1/7/2 neutron or- 
bital) is concentrated in the region where the 1g?3/2 wave 
function is localized. The overlap with this wave func- 
tion is thus the largest and this explains why the filling 
of the neutron 1/7/2 orbital has an important effect on 
the proton 1g?3/2 energy which is strongly lowered. On 
the other hand, when the neutron 2p3/2 orbital is filled 
(from 48 Ca to 52 Ca) the potential changes very little in 
the region where the lcZ3/2 wave function is appreciable. 
This explains why the ld3/2 energy decreases more from 
A=44 to 48 than from A=48 to 52. The energy of the 
2s 1/2 proton state is much less sensitive to the neutron 
shell structure in these isotopes and it decreases rather 
monotonically from A=44 to 52. 



IV. TENSOR FORCE EFFECT 

The tensor force plays certainly a role in the evolution 
of single-particle states. This is discussed, e.g., in the 
framework of the shell model in ref. In a mean 

field approach, the tensor effect originates from the it- 
nucleon and p-nucleon contributions to the Fock terms 




FIG. 6: Potential Phf (see text) (top), square of ld3/2 and 
2sl/2 radial wave functions times r 2 (middle) and their prod- 
uct (bottom) calculated with SkI5 in 44 Ca, 48 Ca and 52 Ca. 



FIG. 7: Difference between the energies of 2sl/2 and ld3/2 
proton states in 40 Ca, 48 Ca, 52 Ca and 70 Ca calculated with 
SLy5 with and without the tensor contribution. 



[3 and it can be introduced phenomenologically 
in the parametrizations of effective interactions built for 
HF models [13, HH, [HJ . Recent progress have been made 
in determining the tensor terms of Skyrme interactions 
[H,H3 by adjusting the single-particle spectra measured 
in N=82 isotones and Z=50 isotopes p5|. 

When this force is included the spin-orbit potential 
presents an additional term depending on the spin den- 
sity J, namely, 



Wn 

U! = ^(Vp + V Pg ) 



aJ q + f3J, 



(5) 



where q stands for neutrons (protons) and q' for protons 
(neutrons), a and (3 consist of a sum of central and tensor 
contributions: a = ac + cxt, /3 = 0c + &t- The central 
contributions depend only on the velocity-dependent part 
of the Skyrme force whereas the tensor contributions are 
gen erated by the tensor component of the Skyrme force 

In, Hi. 

To estimate the effect of the tensor force in our case we 
use the Skyrme force SLy5 which already contains in the 
fitting protocol the terms ac and /3<7-For ar and 0t we 
adopt the values determined in Ref. [23| by comparing 
the Skyrme-HF predictions with the data of Ref. [25| . 
These values are a T = -170 MeV fm 5 and (3 T = 100 
McV fm 5 . We expect that the tensor force favors the 
inversion in 48 Ca (see for instance Fig. 4 of Ref. 0). 
Actually, from 40 Ca to 48 Ca the 1/7/2 neutron orbital is 
filled. The interaction between the proton orbital ld3/2 
and the neutron orbital 1/7/2 is attractive and its effect 
is to lower the energy of l<i3/2, thus favoring the cross- 
ing with 2sl /2. As an illustration we performed SLy5-HF 
calculations for 40 Ca, 48 Ca, 52 Ca and 70 Ca ( 70 Ca is still 
bound within SLy5-HF). We show in Fig 7 the values 
of Ae obtained with and without the tensor contribu- 
tion. As expected, the tensor force increases the slope 



going from 40 Ca to 48 Ca bringing the two states close 
together and improving the agreement with the experi- 
mental data. The improvement is quite strong since the 
value of Ae in 48 Ca is equal to -0.26 McV and -1.41 MeV 
with and without the tensor contribution, respectively. 



V. SUMMARY 

In this article we have analyzed the modification of the 
proton single-particle states 2sl/2 and ld3/2 in Ca iso- 
topes within the non-relativistic Skyrme-HF and the rela- 
tivistic RMF models. Pairing effects have been neglected 
since Ca isotopes are proton closed-shell. We are inter- 
ested in the evolution of proton states and the inclusion 
of neutron pairing does not affect significantly the global 
trend of our results. Both models, HF and RMF, predict 
the same evolution with increasing A for the difference Ae 
of the energies of the two states. This evolution depends 
on the neutron orbitals which are filled, Ae presenting 
maxima and minima corresponding to neutron shell or 
sub-shell closures. In particular, going from 40 Ca to 48 Ca 
the two proton states come closer to each other and they 
can sometimes cross in some models. By performing an 
analysis based on the equivalent potential in the non- 
relativistic Skyrme-HF approach, we have shown that the 
kinetic and spin-orbit contributions present quite a regu- 
lar behavior with increasing A. They both strongly favor 
the inversion of the two states in very neutron-rich nu- 
clei. We have also verified that the contribution which is 
mostly responsible for the maximum of Ae at 48 Ca (and 
leading to an inversion for some models) is the central 
HF potential and, in particular the to and £3 terms. The 
former term favors the crossing of the two states near 
48 Ca whereas the latter acts against it. The net effect is 
that the two states get closer and can cross each other in 
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some models. 

We have finally analyzed the role of the tensor force 
within the SLy5-HF model and found that its contribu- 
tion goes in the same direction as the to term of the HF 
potential, favoring the inversion of the states near 48 Ca. 

Our analysis was restricted to a purely mean field 
picture. Work should be done to include effects be- 
yond mean field. For instance, particle-phonon coupling, 
which has been neglected here, is expected to improve 
the quality of the theoretical predictions in the study of 



single-particle states evolution. 
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